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Q_i' Abstract: We derive analytical expressions for the shape of the invariant mass distribu- 

^ • tions of massless Standard Model endproducts in cascade decays involving massive New 

D . Physics (NP) particles, D — > Cc — > Bbc — > Aabc, where the final NP particle A in the 

' cascade is unobserved and where two of the particles a, b, c may be indistinguishable. 

Knowledge of these expressions can improve the determination of NP parameters at the 
^ ' LHC. The shape formulas are composite, but contain nothing more complicated than log- 

' arithms of simple expressions. We study the effects of cuts, final state radiation and 

detector effects on the distributions through Monte Carlo simulations, using a supersym- 
metric model as an example. We also consider how one can deal with the width of NP 
particles and with combinatorics from the misidentification of final state particles. The 
possible mismeasurements of NP masses through "feet" in the distributions are discussed. 
Finally, we demonstrate how the effects of different spin configurations can be included in 
the distributions. 
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1. Introduction 



While the discovery of a broken TeV-scale supersymmetry ||T|, |2|, ^ ^, ^ at the LHC would 
solve many problems in modern particle physics, it would also raise many questions. For 
example, do the forces unify at some high scale, as the supersymmetric evolution of cou- 
plings imply; what is the mechanism of supersymmetry breaking; does supersymmetry 
provide a viable dark matter candidate? To answer any of these questions it is important 
to have accurate measurements of the supersymmetric partner masses. 

In the Minimal Supersymmetric Standard Model (MSSM), superpartner production 
cross-sections at the LHC are rather large, allowing the discovery of squarks and gluinos 
with masses up to about 2.5 TeV. Nevertheless, if R-parity is conserved, the traditional 
method of measuring masses from resonant peaks in invariant mass distributions cannot 
be easily used for superpartners. The end products of every superpartner decay necessarily 
include the lightest supersymmetric particle (LSP) which is stable and escapes the detector, 
preventing full reconstruction. However, invariant mass distributions constructed from the 
visible particles in a decay chain clearly do depend on the masses of their parents, and it 
should be possible to extract these masses in a systematic way. In particular, the kinematic 
endpoints of these invariant mass distributions (i.e. their minimum or maximum values) 
can be measured and their relation to the unknown superpartner masses can be exploited. 
This method has been widely studied in refs. |P, 0, ^ 11, |l^, p^ . 



This study builds upon the work of ref. [|12[ which investigated the measurement of 
superpartner masses via the endpoint method for the decay^ 

QL X2Q ^R^nq Xiklnq (1-1) 

at the Snowmass mSUGRA benchmark point BPS la [0]. Several problems with the 
endpoint method were uncovered. These include the possibility of multiple minima in 
global least squares fits to the masses and mismeasurements of the endpoints due to the 
occurence of "feet" in the distributions. 

Supersymmetry is not the only New Physics (NP) scenario that can be discovered at 
the LHC and which gives rise to long decay chains. In universal extra dimension models 
(UED) [^], Kaluza-Klein (KK) excitations of Standard Model (SM) particles could have 



experimental signatures very similar to SUSY models. In [16| the possibility of the decay 
chain 

Qi Ziq Lil^q 'jiklnq (1-2) 

was pointed out. Here the subscript 1 denotes the first KK-excitations of the SM particles. 

Indeed the occurence of a similar decay chain in any multi-particle NP model which has a 

good dark matter candidate should not be a surprise. The dark matter candidate should 

be a weakly interacting, neutral and massive particle, thus escaping the detector. This 

particle must be kept stable by some symmetry, in our examples R-parity for SUSY and 

KK-parity for UED, and this symmetry is likely also to conserve the number of NP particles 

in decays, creating a decay cascade through the hierarchy of NP masses. 

^The subscripts n and / on the leptons, representing "near" and "far", denote the first and second lepton 
emitted in the decay. 
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D C B A 

Figure 1: Decay chain. Particles c, b and a are assumed massless. 

In light of this we will maintain a degree of generality in our derivation. We derive 
analytic formulae for the invariant mass distributions of particles in the decay chain 

D ^ Cc^ Bbc^ Aabc, (1.3) 

illustrated in figure ||, where D, C, B and A are massive, satisfying 

mr, > mc > ms > mA, (1-4) 

and the particles a, b, c are taken to be massless. 

With A invisible, we can form four invariant mass combinations of the decay products 
c, b and a: mba, ruch-, mca and nicha- However, both the decay chains (|1.1|) and (|1.2D give 
two opposite sign, same flavour leptons, and Zf, which are in principle indistinguishable. 
This is a general feature of decay chains; given a SM decay product with a particular flavour 
this flavour should be carried on in the decay chain and may result in two same flavour 
decay products. We will deal with this by calculating the shape of the distributions of the 
invariant masses w.c2(iow) and rn^2(h\gh) i which are the invariant masses of c combined with 
either b or a, depending on which gives the lower and higher invariant mass. 

In the following section we derive the shapes of these invariant mass distributions. 
Using the decay chain ( |1.1| ) as an example we consider in section |3| some complicating 
effects present already at the parton level, such as final state radiation, cuts and the 
width of the SUSY particles. In section § we further look at detector induced effects on 
the distributions and in section ^ we discuss a possible application of the shapes of the 
invariant masses to alleviate the problem with "feet" . In section ^ we give a brief summary 
of our results. The inclusion of spin effects is discussed in an appendix. An application of 
the shape formulas to Monte Carlo data for the purposes of improving mass determination 
at the LHC will be presented elsewhere. 

2. Invariant mass distributions 

Our aim here is to obtain expressions for the invariant mass distributions of cascade decays 
such as which can be compared to data to fit NP particle masses. These distribu- 

tions can be easily calculated numerically by Monte Carlo integration, and indeed we use 
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the Monte Carlo event generator PYTHIA [|T^] to numerically check our analytic results. 
However, Monte Carlo integration (via phase space generation) is slow in comparison to 
the evaluation of an analytic expression, which makes it impractical for fitting masses, 
or performing other comparisons, where many distributions for different masses must be 
generated. For this reason we derive analytic expressions. 

In the main part of the paper, we will assume that there are no spin correlations be- 
tween the subsequent decays, and thus treat the decaying particles as if they were scalars 
(see, however, the appendix). Indeed, PYTHIA 6.208 uses this simplification in its gener- 
ation of phase space. This has been shown to be a reasonable simplification for the decay 
chains of physical interest as long as one does not measure the charges of the final state 
particles [|^, |^], which for jets would be very difficult to do. However, it has been sug- 
gested in [|^, ^] that measurements of NP particle spin is possible at the LHC because 
of the asymmetry between produced squarks and antisquarks. The method used in the 
derivation of the invariant mass distributions can fairly easily be adapted to a particular 
spin configuration of the particles in the decay chain and distributions including spin effects 
are given in the appendix for our SUSY example. For simplicity of the analytic expressions, 
we will also ignore particle widths. The effects of widths can be added later by a smearing 
of the distributions. 

Our philosophy for calculating the invariant mass distributions will be to write the 
square of the invariant masses in terms of angular variables in the rest frame of decaying 
particles. In the spin-0 case these have isotropic, flat distributions. After performing a 
transformation from these 'flat' variables to a set of variables which include the invariant 
mass under discussion, an integration over the extra variables can be performed to provide 
the desired distribution. The difficulty of the calculation will be seen not to lie in performing 
the integrals, but in finding the correct integration interval for all possible configurations 
of masses that obey eq. (|1.4| ). Spin effects only complicate the integrand and not the 
endpoints of the integration. 

Distributions formed from two "neighbouring" particles, rricb and ruba, have simple 
triangular shapes when widths and spin are ignored (see e.g. figure 10 and eq. (4.2) of 
fl^), and will not be discussed further here. We will also limit our discussion to cases 
where the intermediate NP particles C and B are on-shell. For a discussion on the shape 
of distributions in the SUSY scenario of eq. (|l.lD , with virtual particles in the decay chain. 



see [21|. 



2.1 The two-particle invariant mass nica 
The invariant mass rrica is given by 
m 



4 = iPc+Pa? = 2i?(^)i^(^) 1 -COS^M , (2.1) 



(B) ( B) (B) 

where Ea and Ec are the energies of a and c respectively and Oca is the angle between 
c and a, all in the rest frame of B, denoted by the superscript (B). Energy and momentum 
conservation in the decays of the B, C and D additionally provide 
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El ' = ^, 2.3) 

^ 2mB ^ ' 

im (771% — m%) rriB 

Ef^ = ^^-^ ^, (2.4) 

where 9^f^ is the angle between c and b in (B). It will be convenient to introduce three 
further quantities, 

t2 \ 1^2 ^2 



(mS^")2 = ^-5 ^ ^, (2.6) 

m 



) 

c 



.2 ^2 ■ 



.max ^2 _ " "^cj K^B " "^a) 7^ 



where and are the maximium possible values of mca and rricb respectively, while 



'^c2{eq) is a possiblc maximum of the mc2(iow) distribution (see sect. 2.3 ) 
From the above, it follows that 

2 ^^ o{B) 

rri n 1 — cos Oca 



Following our philosophy of writing the invariant mass in terms of variables with flat 
distributions we now express cos in terms of cos , the same angle in the rest frame 
of C If C is a scalar this angle will be isotropically distributed. If C is not a scalar, 
the angular distribution will depend on the helicities of b and c. However, summing over 
particles and antiparticles in the final state will cancel these spin correlations and return 
an isotropic distribution. This latter case is applicable here. 
In the rest frame of C we have the familiar result 



ml = (Krf I ^ T"'' 1 • (2-9) 



Using eqs. ( |2.3D and (2^), the same invariant mass in (B) is given by 



(1 - cos 

ml = {m^rf ^ TbS- (2-10) 

(w?(j + ml) — (m^ — m^) cos P 



Equating the two expressions and solving for 1 — cos gives 



1 -COSl/^f, 



\ , m^l + cosg^^^ 
"^ll-cos0^P 



(2.11) 
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and consequently 



(m, 



max\2 
ca ) 



m. 



m, 



C 



^cb 



+ 



1 + COS6' 



cb 



1 - COS ^ " 



(2.12) 



We have now written the invariant mass in terms of quantities which have isotropic 
distributions. Adopting the notation 

_l-cos^ _ l-cosgjf) 

we observe that the differential decay widths for these observables are flat for < (w, v) < 1: 

1 d^To 



Fq dudv 



e{i-u)e{u)e{i-v)9{v), (2.14) 



where 6{x) is the usual step function, and where the subscript "0" is a reminder that 
spin correlations are omitted (however, see the Appendix). We can now make a change of 
variables from {u, v) to {u, rn^) with 

ml, = {mZn\^-au)v, (2.15) 

where, for abbreviation of the equations to follow, we have written 

2 

a = l ^, 0<a<l. (2.16) 



m; 



2 ' 

C 



This gives the differential distribution 

1 a^Fo 5(u, v) 



Tq dudmL 



d{u, m[ 



2 ) 

ca) 



e{i - u)e{u)e{i - v)e{v) 

2 ^ e{u) 



- au) J (m°f^)2(l - au) 



(2.17) 



where for ease of notation we have defined a "top-hat" function 9{x) = 6{x)6{\—x). Finally 
we must now integrate over u to find the distribution of r?T.ca- 

1 aFo r 1 , 

-au 



Ft 

J-oo ^ 



To dml, Fo dudml^ 



"^'^ ^ ' -du 



lo (m-^)2(l-au) 

for < rrica < 'nicct^^ where 

1 



I .^^...^2n ...^ du (2.18) 



^imax = min 1, 

a 



2 



(m; 



max"\2 
ca / 



(2.19) 
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Evaluating the integrals we find 

(2.20) 



1 BTq 



' 1 JTT-B n 

7 775— In — ^ tor < mca < , 

-, max \ 2 



Z In ^" ^ for —m'^^'^ < m < rr?™'^'' 



and zero otherwise. 

From the distribution of the square of an invariant mass m^, it is a trivial task to find 
the distribution of the invariant mass as 

2.2 The two-particle invariant mass mg2(high) 

Having demonstrated the basic method for the simplest non-trivial case, we now turn to the 
more complicated problem of finding the distribution of the observable quantity rnc2{high): 
defined by 

"T'c2{high) = max {rricb, rrica) ■ (2.22) 



Here, m^^ and m^^ are given by eqs. ( ^ ) and ( |2.15 ), respectively. In the notation of 
eq. ( p. 13 ), the required invariant mass squared can be written as 



= max [(mr^)\, {mZ^^ (1 - an) v] . (2.23) 
We now introduce a new variable, 

X = (mS''^)2u - (mJP„^^)2 (1 - au) v, (2.24) 
so that the sign of x picks out which of m-cf, or rrica is larger. Then 

™c2(high) = 0{x){m2n'u + ei-x){mZn' (1 - au) v. 
The new variables can be inverted to give 



"^c2{high) + ^(-^)^ _ "^c2(high) - ^ (2;) X 



(2.25) 



Since the double-differential width with respect to u and v is flat [see eq. ( |2.14 )], we may 
write 



1 d^To 
To dxdml^^^.^^^ 



d{u, v) 



d{x, m? 



c2{high) 

e{u)e{v 



e{u)e{v) 

' (2.26) 



(m^a^)2(m;5^^)2(l - auY 
and integrate over x to give the desired distribution: 

1 dVo f^^ ^ 1 

To 5m2^(high) ^""^^ ^'^^ )2(m--)^(l - an+) 
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+ 



J — oc 



(n_)0(t;_) 



-dx. (2.27) 



(mmax)2(-^max)2(l - au_ 

In the above we have written u± and v± for u and in the case of positive/negative x, i.e., 

+ X m 



m 



c2(high) 



c2(high) 



(m--)2 1 - a 



c2(high) 



(^max)2 ' 
,2 



"^c2(high) ^ 



(m--)2 1 - a 



j£2{high) 



(2.28) 
. (2.29) 



While the integrals themselves are trivial, we must take special care with the integration 
endpoints. The step-functions restrict < U-t < 1 and < U-t < 1, which in turn give 
restrictions on x and/or ?Tic2{high): 



Q(u^) + 

Q(uJ) / 
Q(yJ) / 



< mc2{high) < "^^h""^ 



+ 1 



"cb 



m 



c2(high) 



("^m^'') < ^ < ''^c2{high) ' 



-m 



<x < {m^^^'f 



c2{high) 



m 



c2(high)' 



X < 



[m 



max\2 
cb ) 



(mg^")2 



max") 2 
ca ) 



+ 1 



m 



c2(high)- 



(2.30) 
(2.31) 
(2.32) 
(2.33) 



The first two inequalities constrain the integration over positive values of x, while the last 
two constrain the integration over negative values of x. Notice that V- only provides one 
inequality since the condition 9{v-) ^ holds when (2.32) holds. Also 0[vj^) ^ yields an 



upper bound on x because the denominator of f+ is always positive when ( 2.30 ) holds. 



Which of these inequalities provides the strongest bound and thereby the endpoint 
of the integration is highly dependent on the mass hierarchy between particles A, B and 
C . In particular, the various bounds on x coincide at four (non-trivial) distinct values of 
mc2(high)) namely m™*^^, m^^, ^^f^^a^ "^S'(eq)' ^"^^ relative size of these 

four quantities which is important. With this in mind we define three different regions, 
exhausting all possible hierarchies. Writing 



m . 



mi 



m 



B 



m: 



c 



Ra = — f-, Rb = — f-, Rc = — T' 



(2.34) 



D 



{Rc is defined for later reference), the different mass hierarchies can be divided into: 



Region 1 
Region 2 
Region 3 



2-Ri 



< Rb < I, 



Ra < Rb < 
< Rb < Ra- 



The three regions are shown in figure over the space of Ra and Rb- 
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2.2.1 Region 1: < Rb < 1 

In Region 1, we have 



max ^ ™max , max ^ _„max ocN 



The integration over positive values of x is restricted by eq. ( 2.31 ), and is non-zero only 



when ?Tij;2{high) satisfies eqs. ( 2.3C1| ). If the lower bound on x is greater than zero, it will 



provide the lower limit of the integration. However, this is only the case when 

rnc2(high) > "^S(eq) (2-36) 



which can never be true for Region 1 due to eq. ( p.30| ), and there is only one permitted 
range: 



< X < ml^^high) for < m,2(high) < "^r^"- (2-37) 



For the integral over negative values, x is restricted by eqs. (|]3|) & (|^. While the 
lower bound is unambiguous there are three possible upper bounds (including the original 
X < from the integral) and we must divide the solution into three possible cases depending 
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on which upper bound is dominant: 



-"^c2{high) <X<0 



for < 



"^c2(high) < ITT-cb ' 



-m 



c2{high) 
2 



m 



(^max)2 



""^c2(high) < X < 



c2{high) 



for m^^- < me2(high) < ^m^''", 



+ 1 



c2(high) 



for ^mZ^'' < mc2{high) < m, 



max 
ca 



(2.38) 



Performing the two integrals in the appropriate cases and adding them together we 
find the full expression for the "T'c2(high) distribution: 



1 dTp 

To 5"lc2(high) 



(^max)2 



In 



am 



+ 



c2(high) 



■In- 



m; 



C 



I ("^?a'^")'a "^c2(high) 



and zero otherwise. 



for < mc2(high) < m 
for m™''^ < mc2(high) < 



max 
c6 ' 



max 



for ^"^m " < "ic2(high) < 



m" 



(2.39) 



2.2.2 Region 2: Ra < Rb < 2^ 

For Region 2 we have 



mc 



max ^ ^max ^ ™max ^ ^max 
"T-ca < "^c2(cq) ^ ^cb < 



(2.40) 



As for Region 1, the integration over positive x is restricted by eq. ( |2.31 ). However, unlike 



Region 1, the lower limit of eq. (2.31) can now be larger than zero, providing a new lower 
limit of integration. This happens when eq. ( p. 36 ) is satisfied, and provides different limits 
of integration for two distinct cases: 



< X < m^2(high) for < "^c2{high) < "i;;2(eq), 



a(m™'''')2 



+ 1 



"^c2{high) 



[m 



max\2 



r <x< ml^. 



(high) 



for 



c^cq) < "^c2(high) < 



(2.41) 



For the integration over negative x, the step-functions again require that eq. (|l3^ ) 



and eq. ( |2.33| ) hold, if the integral is to be non-zero. However, since is larger than 



ms.™max 



in Region 2, the upper bound from eq. ( 2.32| ) is never dominant. Therefore 
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Region 2 gives only two different sets of integration limits: 



-m 



^ < X < 



c2(high) 



""^c2(high) < ^ < 



for < mc2(high) < rn^2'' 



(eq)' 



(^max)2 



+ 1 



■c2(high) 



for m" 



c2(eq) < "^c2(high) < "^ca 



(2.42) 



For Region 2 the full distribution is then given by 



1 arp 

To 5"ic2(high) 



1 
1 



(^max)2 



In 



am 



+ 



c2(high) 



In 



c2{high) 



+ 



m 



c2(high) 
('™max\2 



m 



max'\2 
c6 I 



and zero otherwise. 



for < mc2(high) < 



for < 



m 



c2(eq) ' 



c2(eq) < "^c2(high) < "1^^ , 



for m^''^ < mc2(high) < 



(2.43) 



2.2.3 Region 3: ^<Rb <Ra 



In Region 3 we have 



— — < "ic2(cq) < "^ca < "T-cb • 



(2.44) 



For the individual integrations over positive and negative x, Region 3 is identical to Re- 
gion 2, and we again obtain eqs. ( |2.41| ) & ( 2.42| ) for the positive and negative integration 
limits respectively. However, since m^^^ is now larger then rn^^, contrary to Region 2, 
the sum of the two contributions will be different. We find: 



1 dVo 

To 5?"c2(high) 



{m^^)'^a 



(^max)2 



In 



am 



+ 



c2(high) 



(^max)2 _ am22(i^jgj^) (mJ?^^-)2 - am^^^^.^^^ j 



for < mc2(high) < 



m 



c2(eq) ' 



1 



(m--)2a 
1 



"^c2(high) ' i^c^y 



for "^S'fcq) < "^c2(high) < 



m' 



and zero otherwise. 



for mfj''' < ?nc2(high) < "^S"'''' 

(2.45) 
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2.3 The two-particle invariant mass m(,2(iow) 
The invariant mass m^2{\ow) is given by 

"T'c2(low) = min {rricb, ruca) ■ 

For calculation of the differential decay rate we can adopt a method very similar 
of ?7T,c2(high) above. Using u and v as defined in eq. ( 2.13D we can write ?7t,c2(1ow) as 

"ic2(iow) = [imT"?u, {ml 
Keeping x as defined in eq. ( |2.24| ), we write 



,max\2 



1 — au) v\ 



max\2 



These can be inverted to give 
u 



"^c2(low) + ^(2^)^ 



m 



c2(low) 



{—x) X 



^^max)2 



"^c2(low)+^(^')^ 



(2.46) 

to that 

(2.47) 
(2.48) 

(2.49) 



and the differential distribution again has the form ( p.27| ), but now with 



m 



[m 



c2{low) 

max\2 ' 
cb ) 



"^c2{low) + ^ 



(^max^)2 



c2{low) 



X 



m 



c2{low) 



(m--)2 1 



"^c2(low) 



(m--)2 



1 _ „ e2(low) 



(2.50) 
(2.51) 



These are actually the same definitions as eqs. ( |2.28 ) and ( 2.29| ) only now written in terms 
of "ic2{iow) rather than mc2(high), using 

(2.52) 



II 2 2 

Pl = "^c2(high) ~ ™'c2(low)' 



as follows from ( p. 23] ) , ( p.24| ) and ( |2.47| ) . The step-functions of the integrand then give the 
following restrictions on x and/or mc2(iow)- 



0{n+) + 
Q{v^) + 

d{uJ) ^ 
e{v-) ^ 



-m 



c2(low) < ^ < ("^"D - "^c2(low) 



(^max)2 



X < 



(^max)2 



aim 



max "12 
ca ) 



+ 1 



m 



c2(low)' 



< mc2{iow) < m_ 



cb ' 



a{m 



max'\2 



(^max)2 



+ 1 



m 



c2(low) 



(m, 



max\2 



< X < m 



c2(low) • 



(2.53) 
(2.54) 
(2.55) 
(2.56) 



The parallels with the ?7ic2(high) case are obvious. The step-functions provide the same 
constraints for mc2(iow) as they did for mc2(high) ) but constraints which were previously for 
negative x are now for positive x and vice versa. Therefore we must be careful with the 
X = boundary: we will have the same regions of applicability as shown in figure ^ but the 
inequalities giving the integration limits for each region will be different from the ?7ic2(high) 
case. 
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2.3.1 Region 1: < Rb < I 

In Region 1, eq. ( |2.35| ) liolds as before. 



For positive x, we have constraints given by eqs. ( |2.53| ) & (2.54). However, in Region 1 



the upper bound on x provided by eq. ( p.53|) is always more restrictive, and insists that 
mc2(iow) be smaher than for a nonzero result. Finally, since "^c2(low) necessarily 

positive, the lower bound on x is trivially zero. Therefore, the integration limits for Region 1 
with positive x only has one case: 



< X < {mfr f - m^^d^^) for < m,2(iow) < "^S""- (2-57) 



For negative x, the integration limits also only have one case since eq. ( |2.55| ) restricts 
"ic2{iow) to be below all the characteristic masses in eq. ( 2.35 ). This gives 



-(m, 



maxN2 



+ 



m 



■c2{low) 



< X < for < mc2(iow) < fK^b 



(2.58) 



It is then rather trivial to calculate the full integral for Region 1: 



1 



Fn dm?r,n ^ 
" c2(low) 



max "12, 



In 



(m, 



max\2 
cb ) 



am 



c2(iow) a[m 



max^2 



am 



.(^max)2 



+ 



c2(low) 



max\2 
cb ) 



^j^max^2 



am 



c2(low) 

(2.59) 



for < mc2(iow) < fnf^^, and zero otherwise. 



2.3.2 Regions 2 and 3: < < 

The inequalities given by the step-functions together with the positivity or negativity con- 
straints on X, are independent of the boundary between Regions 2 and 3. In other words, 
which of m^^ and m^^^ is larger is irrelevant when both are larger than nT-^^^y Therefore 
the analytic form of the distribution will be the same in Region 2 as it is in Region 3 and 
we do not need to treat them separately. The hierarchy of characteristic masses is then: 



mc 



max ^ ™max 
^ca ^ "^c2(oq) 



< min(m™'''',mJ5'''') 



(2.60) 



For positive x, the issue is which of the upper bounds of eqs. ( 2.53 ) & ( p. 54 ) is more 
restrictive. For higher values of w.c2(iow) it is eq. ( p. 54 ) which is more restrictive, with the 
transition being at ^^m^^ as already intimated by eq. ( 2.38| ), and we must also ensure 
that this upper bound on x is positive. Together, these considerations give the integration 
limits: 



< X < {m"}^^)"^ - m^. 



< X < 



'cb 



c2{low) 

■ (mmax)2 



max 



+ 1 



"^c2(low) 



for < mc2(iow) < 



for ^m-- < mc2(iow) < mf^^^y 



(2.61) 
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For negative x, only the (lower) constraint of eq. ( 2.56 ) is interesting. The insistence 



that this lower bound is negative requires that m^2(iow) is smaller than m-J^^'^gq), which is 



more restrictive than eq. ( |2.55 ) for these regions. The integration limits for negative x are 
therefore again rather simple: 



(mr^)2 + 



max\2 



(-^max)2 



+ 1 



"^c2(low) <X<0 



for 



< m,2(low) < m;^2{eq)- 



(2.62) 



The full result for Regions 2 and 3 is: 
1 dVo 

,max\2 



(m--)2a 



1 



111 Q ~|~ 



am 



c2(low) 



B /'™max\2 



/'TOmax\2 /'™max\2 ,im2 
l''*c6 'I V"*cfe ^ ""*c2(low) 

"^c2{low) < ^ 



for < m.on.w^ < ^m!!i'^" 



max\2 



In- 



(m; 



am 



c2(low) 



(^max)2 (TOmax)2 



am 



("^r")'"^c2(low) 



+ 



c2(low) 



am 



for ^mf^"" < mc2(iow) < "i""''' 



c2(eq)' 



c2(low) 



(2.63) 



and zero otherwise. 



2.4 The three-particle invariant mass mcba 
The quantity mcba is defined by 

^Iba = (Pc+Pb+Paf- 



(2.64) 



Making use of energy and momentum conservation, we can write this in terms of particle 
masses and two angles: 



"icfea = m^, - m^ - y 



m^ + m' 
rriB 



+ Jy^ + 2ymB + m| - mj-, 



2 ^2 / "^B "^A 



(1 - 2u;) , 



(2.65) 



where u is given by eq. ( 2.13 ), w is similarly defined in terms of the angle between A and 
D in the rest frame of B, 



1 — cos 



w = 



(2.66) 
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and the quantity y is defined by 



2ml 2mBml \ 2mB ' ^ ^ ^ ^ 



B 



It is straightforward to solve for u and w, 

{ml — ml^ ml — 2mlmBy 
{ml - ml) {ml - ml) 



mD-m%-y[ 1 + Jy^ + 2ymB + mj^ - m^ ( '"-^^"^ ) - m 



'■cba 



w = ^ ; ^ ' " ' . (2.68) 

2 J y2 + 2ymB + m| - m^, ( ^^^^ 



Under the assumption that i? is a scalar particle, its decay is istotropic and the doubly- 
differential width in terms of these two variables u and w must be flat. Therefore we can 
write 



1 







To dydm^ 



cba 



d{u, w) 



d{y, m^ 



cba J 



1 92r, 







To dudw 



9{u)9{w) 



[m 



D 



ml) {ml - ml) {ml - m\) ^y2 + 2ymB + ml - ml 
and integrate over y to acquire the desired distribution, 



(2.69) 



1 dVn 



Tq dm 



cba 



{ml 



6{u)e{w) 



a J-oo 



y2 + 2ymB + ml — ml 



--dy. 



(2.70) 



At this point it is useful to make some definitions, in order to help keep expressions 
compact. Let us first define the integral: 



L{ai,a2) 



""^ \l + '^yn^B + ml — ml 



--dy 



In 



a2 + mB + a/«2 + 2a2mB + m^ - mj, 



(2.71) 



ai + mB + a/ «i + 2aimB + ml — ml 



which is the only integral that we will need. The mass values which appear as kinematical 
endpoints of mcba are |1C, 12 1: 



m\ 



[m 



D 



d 



m 



(2.72) 



c 
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^2 _ - ^b) [mh^B - mlm\j 
'"2 ~ _2 „„2 ' 





2 

C 


-ml) 








1 1 








2 






D 


-ml) 


{ml - 



) 



ml = ^ ^ ^, 2.74 

ml 

ml = {mo — niAf ■ (2-75) 
The mass m4 is always the largest, as is seen by 

2 2/2 \2 ; 2 

771-4 — m-^ = [mQ — mAmu) /'mQ, 

m\ — m\ = (m^m/) — m^m^)^ /(m^m^), 

m\ — m\ = [ml — m Am d)"^ /ml- (2.76) 

For the others, the relative order will depend on the values of mA, ruB, mc, and mB- 

m\-m\= [ml - m^) (m|m|, - m\;m\.) /{mime), 
m\-ml= {ml - ml) {m\ml - m|m^) /{mlml;), 

ml-ml= {ml - mc) {m\mc - mlml) /{mime). (2.77) 



As for the previous distributions, it is important to determine the upper and lower 
integration bounds. The step-functions for u provide rather simple constraints, 

yi = — ^ ^ ' <y< — ^ = y2, 2.78 

2mB 2mB 
where we have defined the quantities 7/1^2 for later convenience. The constraints from the 
step- functions for w are somewhat more complicted, 



where 



y3<y< yi, (2.79) 

^^'^ ^ 4mlm\ { ^"^^ ~^a)^ ~ [^cba - {ml - ml)] {ml + m\) 

± {ml - ml) ^J {ml -m\- mlj^ - 4m\ml,^ | . (2.80) 



In order that 7/4^3 be real, we must insist that m^ba < with 777,4 defined by eq. (2.75). 

We now need to compare these constraints and determine the relative ordering of 7/1, 
7/2, ys and 7/4, to see which will form the endpoints of the y integration. There are four 
possible cases where the constraints overlap, as illustrated in figure |3|. 

7/1 vs. 7/4: In order to have a physical solution, we require that yi < 7/4. This con- 
straint is manifest differently for different mass hierarchies. Comparing the expressions 
from eqs. ( p.78| ) & ( p. 80 ) we find the requirements. 



for ml < mATUD (^ mAml < mlmo), yi < 2/4 if rucba < mi, 

for ml > mATTiB & mAml > mlmo, yi < 2/4 if mcha < ^2, (2-81) 

for 777^ > mArriB & mAml < mlmo, yi < 2/4 if mcha < m^, 
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where mi, m2 and 7714 are given by eqs. (|]72D, (|]7|) & ( ^) respectively. 



(a) 

(^) 
(7) 
(•5) 



yi 



y2 



ys 
yi 



ys 



yi 



ys 



y2 



y4 



ys y4 
yi y2 



y4 
y2 



y4 



y2 



y4 



yi 



y2 



y4 



Figure 3: Schematic representation of integra- 
tion ranges. The cases {a)-{6) give non-zero 
contributions. 



2/2 vs. 7/3: Similarly, to provide an overlap of 
the integration regions we must insist that 
2/2 > 2/3- Once again, which values of nicba 
are required is dependent on the mass hierar- 
chy: 

for m| < uiAmD, 2/3 < 2/2 if rricba < m^, 
for m| > uiAmD, 2/3 < 2/2 if rricba < m^, 

(2.82) 

where ms is given by eq. ( 2.74| ). 



7/2 vs. 7/4: Although one does not require any 
particular relation between 7/2 and 7/4 in order 
to have a non-zero result, only the smallest of 
7/2 and 7/4 will provide the upper bound on the 
integration, so it is important to verify which 
mass regimes lead to which dominant upper 
bound. We find: 



for < mArriD, 
for TTT,^ > mAin£), 



2/2 < 2/4 if rricba < ?"3, 

2/2 < 2/4 if rricba < ?"4- 
(2.83) 



7/1 vs. 7/3: Finally, the lower bound of the integration is governed by whichever of yi and 
2/3 



IS 



larger. We find: 



for ml < TJiAniD (^ m^m^, < m|m,,), yi > ^3 if ^2 < m^a < m^, 



for > mArriD 
for > mAiTiD 



ruAml > mlrriD, 2/1 > 2/3 if mi < m^, < m,, 



2/1 



(2.84) 



> 7/3 if min(mi,m2) < rricba < max(mi,m2). 



As for the simpler two-particle invariant masses, the mass-dependent relations ob- 
tained above result in different distributions in different regions of mass-space. We must 
therefore specify the mass hierarchy of these regions before continuing. Given the ordering 
of eq. (|1.4|), an unambiguous division^ is: 



Region 1 
Region 2 
Region 3 
Region 4 



mc < niAmD, 



22 2 
nT'AmQ < m^mD, < niAmD 



ml > mATTiD, rnATTil > m^mo, rn\ < ^A^d 



2 ^ 2^2 
mQ > mAmo, ^A^c < itt-b''^d. 



> mAinD 



ml>mAmn, mAml<mlmo, m| < m^mz. 



(2.85) 
(2.86) 
(2.87) 
(2.88) 



"Note that the 'regions' defined here are diflferent from those of sects. 2.2 and 2.3. The present ones are 



three-dimensional volumes, whereas those of sects. 2.i and 2.3 are two-dimensional areas 
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Figure 4: Allowed values of Ra, Rb and Rc for the four main regions of rricba in terms of 
volumes inside the unit cube. The allowed values for Region 1 lie in the volume below Surface I: 
< Rc < RaRb < 1- For Region 2 the allowed volume is that above Surface II: < ^ < Rc < 1 
(not including the flat region where Rb > Ra)- Likewise the allowed volume for Region 3 is that 
above Surface III: < < Rc < 1. Region 4 lies in the volume above Surface I and below 
Surfaces II and III. (See figure [l^ for projections onto the i?yi-i?B-plane). 



This division, in terms of the mass ratios of eq. ( 2.34 ), is illustrated in figure ^, and has 
the property that each region has a different endpoint for the invariant mass distribution. 
The different orderings of mi, m2 and ms, as given by eq. ( |2.77 ), will lead to further 
subdivisions of the above regions. 

2.4.1 Region 1: < Rc < RaRb < 1 

The first inequality (underlined) of eq. ( 2.85| ) implies the other two hold for our mass 



hierarchy. Also, as seen from (|]7D and ([ITTI) , since < mAniD < msmD and rn^rn^ < 



uiq < m^m|), we have m2 < mi < and m^ < mi < m4 respectively, but the relative 
size of m2 and ma is undetermined. Thus, we must divide the region into subregions. 



Region (1,1): If m^ > mAmc then m2 < m^ < mi (see eq. ( 2.77D ), and we find the 
invariant mass distribution is, in terms of the function L defined in eq. ( |2.71j ). 



1 dVo ^ 1 
^0 9^2. (m^'^^)2a 



L{y3,y2) for 0<mcba<'m2, 
L{yi,y2) for m2 < m^ba < ms, 
^L{yi,yA) for m3 < mcba < mi, 



(2i 



and zero otherwise. 

Region (1,2): Alternatively, if < mAmc then m^ < m2 < mi and the invariant mass 
distribution is: 



1 dT 







To 9m2j,^ 



{m^^^ya 



Hy3,y2) for 0<mcba<m3, 
L{yz,yA) for m3<mcba<m2, 
^L{yi,yA) for m2 < mcba < mi, 



(2.90) 
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and zero otherwise. 



2.4.2 Region 2: < ^ < Rc < 1 

In this region, the second inequahty (under Uned) of eq. ( p.86 ) impUes the other two. Also, 
since niBmi) < and < j^m'^ < mAfnc we have mi < m2 and < ?n2, but the 
relative magnitude of mi and ms is undetermined and we must divide our region in two. 



Region (2,1): If mAmi) < msmc then mi < m^ < m2 (see eq. ( p.77D ), and we find the 
invariant mass distribution is: 



1 dT 







To dm^ 



cba 



{m^^ya 



L{y3,y2) for 0<mcba<mi, 
L{yi,y2) for mi < mcba < ^3, 
^L(yi,y4) for m^, < mota < m2, 



and zero otherwise. 



(2.91) 



Region (2,2): If mAm£, > msmc then 7713 < mi < m2, and the distribution is: 
1 dVo 



To dml^ 



cba 



1 

{m^^^Ya 



L{y3,y2) for 0<mcba<ms, 
L{y3,y4.) for m-i < mcba < mi, 
L{yi,y4) for mi < mcba < m2, 



and zero otherwise. 



(2.92) 



2.4.3 Region 3: < ^ < Rc < I 

In this region, it is the third inequality (underlined) which leads to the other two. The 
inequalities mAmc < mAmi) < and mAm,D < msmc tell us that m2 < m^ and 
mi < 7723. Although we do not know the order of mi and m2, the endpoints of the 
integration are independent of which one is larger and their relative magnitude only affects 
the region of applicability of the different functions (unlike in Regions 1 and 2 where the 
integration endpoints also changed). For the entirety of Region 3 we find: 



1 dTn 



1 



max \ 2 ^ 



L{y-i,y2) 
L{yi,y2) 

^L{ys,y2) 



for < mcba < min(?TT,i, 777,2), 

for min(77ii, 7772) < 777c;,a < max(m.i, m.2), (2.93) 

for max(777i, 7772) < 777cfea < "73, 



and zero otherwise. Note that 7771 < 7772 when msmu < m^, from eq. ( 2.77 ) 



2.4.4 Region 4: < RaRb < Rc and Rc < ^ < 1 and i?c < ^ < 1 

Region 4 encompasses all the other allowed areas of mass-space. In this region we only 
know, a priori, that 7771, 7772 and 7773 are smaller than 7774 (as for all regions) but not their 



Ra 
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relative sizes. However, as for Region 3, the ordering of mi and m2 does not change the 
endpoints of the integration, but only changes the region of applicability. Therefore, we 
have three possible subregions, corresonding to the relative sizes of ms and the maximum 
and minimum of mi and m2. 

Region (4,1): If m^m/j > m^mc and m^ < m^m^ then < min(mi,m2), and the 
invariant mass distribution is: 



1 dVn 



and zero otherwise. 



L{yz,y2) for < mcfea < ms, 

L{y3,y4:) for ma < rricba < min(mi,m2), 

Liyi^yi) for min(mi, m2) < mcfea < niax(mi, m2), 

^L{yz,y4) for max(mi, m2) < mcha < "14, 



(2.94) 



Region (4,2): If {m^m^) > m^, uiAmD > mBmc and m^ > m^mc} or {niBrriD < rriQ, 
mArriD < mBmc and m^ < m^mc} then min(mi,m2) < m^ < max(mi,m2), and the 
invariant mass distribution is: 



1 dTp 
To dml, (m 



1 



max\2 , 
ca ) 



'L{y3,y2) for < mc^a < min(mi, m2), 

L{yi,y2) for min(mi,m2) < m^ba < 1^3, 

L{yi,yi) for ma < mcfea < max(mi,m2), 

,L{y3,y4,) for max(mi, m2) < mcba < "14, 



(2.95) 



and zero otherwise. 



Region (4,3): Finally, if myim/5 < mBmc and m^ > mAmc then max(mi,m2) < ma, and 
the invariant mass distribution is: 

'L{y3,y2 

1 dVo 1 



^0 9m2 (mS^^)2o 



and zero otherwise. 



for < m^ba < min(mi,m2), 

L{yi,y2) for min(mi, m2) < m^fea < rriax(mi, m2) 

L{y3,y2) for max(mi,m2) < mcfea < "la, 

.L{y3,yi) for ma < mcba < ^,4, 



(2.96) 



3. Parton level 

In order to explore how the derived expressions can be applied to real data, we will compare 
with SUSY cascade decays in Monte Carlo (MC) events generated for the mSUGRA Snow- 
mass benchmark point SPS la |14]. SPS la has the GUT scale parameters mo = 100 GeV, 
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mi/2 = 250 GeV, = -100 GeV, tan/? = 10 and /i > 0. The TeV-scale SUSY mass 
spectrum is calculated from these parameters, together with a top mass of mt = 175 GeV, 
by ISAJET 7.58 ||22| (for numerical values, see, e.g. [^]). We have generated a number of 
events equivalent to 50 fb"^ at the LHC, by running PYTHIA 6.208 with CTEQ 5L 
parton distribution functions |23|. 

In particular, we focus on the decay chain ( |1.1D , which was previously investigated in 
some detail in IQ. The comparison is done for the invariant masses mg/(iow)) "^-^/{high) ) '^qif 
and mqii.^ For the quark q we consider only the up quark and the corresponding squark. 
The complications introduced with more than one flavour, and thus multiple, overlaying 
distributions, will be faced in the next section. 

We find the analytic expressions for the distributions by going from the distribution 
of the square of the invariant mass to the distribution of the invariant mass by eq. ( 2.21] ) 
and substituting A = Xi, B = Ir, C = X2 ^^'^ D = ul in our expressions. 

However, this simple picture is not what will be seen in an actual experiment. The 
shapes will be distorted by the width of the SUSY particles, cuts applied to remove back- 
grounds, by Final State Radiation (FSR), and by detector effects. In order to use the 
shapes for extracting NP particle masses from real data it is important to understand 
these effects, and the limits they place on the use of the shapes. Here we consider the three 
main effects apparent at parton level, width effects, the bias introduced by cuts and a shift 
in invariant mass from FSR. 

3.1 Width 

Since we have not taken into account the width of the decaying particles we addition- 
ally smear the analytic expression with a (truncated) Breit-Wigner function. If g{m) is 
the original distribution we plot in the following figures the function f{m) given by the 
convolution 

Ts rm+iOTs g{m') , 

•^^"'^ = 4arctL20y„_ior. (m' - m)^ + (r,/2)2 ^^'^^ 
where F^ is a parameter of the smearing, determined by a fit to the MC data. We normalize 
the analytic expression to the number of MC events. 

The smeared analytic curves (for the SPS la parameters) are compared with the parton 
level MC distributions, with no cuts or Final State Radiation, in figure ^. The agreement is 
very good, being at the level expected from statistical fluctuations in the MC, and provides 
a verification of our analytic expressions. The smearing is responisble for a rounding of the 
sharp peaks in the m.gi(iow) aiid racha distributions, and greatly improves agreement near 
the upper end of the distributions. 

3.2 Cuts 

We apply the cuts used in |]l^] to isolate the decay chain from Standard Model background, 
• the three hardest jets have pj^* > 150, 100, 50 GeV, 



^We include the distribution of vriqij for comparison with our analytical expression, even though it will 
be difficult to separate the near and far leptons experimentally. 
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Figure 5: Invariant mass distributions at parton level with no cuts and no FSR. MC events (black) 
are shown with errors. The lines (red) are the analytic expressions with nominal masses. The dashed 
lines show the ranges of validity of the different function pieces. 

• ET,miss > max(100 GeV, 0.2 Mes), where Mes = -ET.miss + Y1^=iI^t,v 

• the two hardest leptons have p^^^ > 20, 10 GeV. 

Figure ^ compares the analytic curves with the parton level MC events after the application 
of these cuts. From a visual inspection we see a fairly good agreement between the MC 
events and the analytic functions, except that the number of events at low invariant mass is 
reduced in the MC distribution compared to the analytic curve. This is most pronounced 
for m,i(io„). 

In order to determine which of the cuts is introducing this discrepancy, we show in 
figure H the average value of the cut variables -BT.miss — max(100, 0.2Mefj), Pj'\, Pj^2 ^^'^ 
p^j'\, in each bin of the mg^pow) invariant mass distribution versus the invariant mass in 
that bin. The two jet cuts omitted have a very similar behaviour to that on the hardest 
jet. 

For the cut on missing energy there is a large gap between the signals's missing energy 
and the cut value (except for the two end bins, where low statistics lead to large errors) 
and therefore the effects on the shape of the invariant mass distribution will be small. On 
the other hand, because of the small difference between the cut value and the average value 
at low invariant masses, the cut on the transverse momentum of the second hardest lepton, 
p^'^, will clearly introduce a bias by removing more events at low invariant masses. From 
figure 1^ we see that we should be safe in trusting the analytic shape of the distribution 
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Figure 6: Invariant mass distributions at parton level with cuts as described in the text, but no 
FSR. MC events (black) are shown with errors. The lines (red) are the analytic expressions with 
nominal masses. The dashed lines show the ranges of validity of the different function pieces. 

after the cut on p^^, down to around niqKiovf) = 100 — 150 GeV where the distribution 
flattens out due to the cut imposed. The cut on the transverse momentum of the hardest 
lepton, Pj^i, will likewise have some effect on the distribution, but only at very low invariant 
masses. For the cut on transverse momentum of the hardest jet, and similarly for the other 
two jet cuts, we again find that the average values do not lie close enough to the cut value, 
for most of the invariant mass bins, to cause a problem. However, even if this had not been 
the case, since the average value of as a function of is to a good approximation 

flat, it would have cut equally at all invariant mass values, and so introduced no bias. 

Studying the average of a cut variable over the range of possible invariant mass values 
enables us to give an estimate of where it is safe to fit the distributions given some cut 
value, and could also be used to optimize the cuts so that they introduce less or no bias 
for certain ranges of invariant mass. 

3.3 Final state radiation 

In figure ^ we include FSR in the MC data and we observe a slight shift towards lower 
invariant masses for all four distributions as compared to the analytic shape. This can be 
compared with figure |6| which included no FSR. 

Ignoring statistical fluctuations between the two different event samples this shows 
that FSR conserves the shape of the distributions, but lowers the invariant masses slightly, 
indicating that it should be possible to use the shape of the distributions in determining 
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Figure 7: Average value of the cut variables i?T,miss — niax(100, 0.2Moff), p^^P'i, p!^*^ andp^^*j^, versus 
invariant mass for rrigKi^^-^ (black). We impose the cuts described in the text on the invariant mass 
distribution before computing averages, but include no FSR. The horizontal lines (red) show the 
cuts described in the text. 



SUSY masses, but that the precision will be limited by that of the jet energy scale. 



4. Detector effects 

In section ^ we showed that the shape of the invariant mass distributions from the decay 
chain in eq. ( |1.1[) to a large extent survived parton level effects and cuts to remove SM 
background. We demonstrated a method of showing whether and where cuts could deform 
the distribution, making our shape predictions from the kinematics of the decay chain un- 
reliable. What remains to be discussed are the effects of the detector in a given experiment 
and the combinatorical complications introduced by trying to pick jets that correspond to 
the quarks in the decay chain. 

For this purpose we use AcerDET 1.0 [^], a generic fast detector simulation for the 
LHC, similar in structure to the ATLFAST |25] MC simulation of the ATLAS detector. 



AcerDET expresses identification and isolation of leptons and jets in terms of detector 
coordinates by azimuthal angle (j), pseudo-rapidity r] and cone size AR = \J (A0)2 + (Ar/)^. 
We identify a lepton if > 5(6) GeV and |?7j < 2.5 for electrons (muons). A lepton is 
isolated if it is a distance R > 0.4 from other leptons and jets and the transverse energy 
deposited in a cone AR = 0.2 around the lepton is less than 10 GeV. Jets are reconstructed 
by a cone-based algoritm from clusters and are accepted if the jet has > 15 GeV in 
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Figure 8: Invariant mass distributions at parton level with FSR included, but no cuts. MC events 
are shown with errors (black). The lines (red) are the analytic expression with nominal masses. 
The dashed lines show the ranges of validity of the different function pieces. 



a cone AR = 0.4. The jets are recalibrated using a flavour independent parametrization, 
optimized to give a proper scale for the dijet decay of a light (100 GeV) Higgs particle. 

As in the parton level discussion we use the mSUGRA model point SPS la, a number 
of events equivalent to 50 fb^^ of integrated luminosity generated by PYTHIA 6.208 |17|, 



and the cuts of subsect. 3.2. In addition to the cuts introduced earlier we also cut on 



6-tagged jets to remove events with a 6-squark in the decay chain, which will, due to a 
smaller squark mass, have a different distribution. We assume 50% 6-tagging efficiency with 
rejection factors 100 and 10 on jets from gluons/three lightest jets and c-jets respectivly. 
While this is on the conservative side it removes a majority of events with 6-squarks in 
the decay chain. However some events remain that will subtly change the distribution. 
How large a problem this will be will depend on the rate of b production and the efficiency 
achieved for 5-tagging. We have the same issue for the u and d-squarks, since they in 
general have different masses, but here we cannot tag the jets. One possible solution to 
this problem is to fit the experimental distributions with the weighted sum of two functions 
with different squark mass values. This effect is of course smaller than for the 6-squark, as 
the mass difference of these squarks is relatively small in our scenario. 

Even if we could remove effectively all events with 5-jets another potentially large 
problem remains. We have no sure way of knowing which jet corresponds to the quark in 
our decay chain. Since the other jets can stem from the decay of the other SUSY particle 
produced in the hard process or from the underlying event we cannot rely on "high-low" 
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Figure 9: Invariant mass distributions after detector effects and consistency cuts shown with 
error bars (black). Also shown, are the parton level distributions (blue histograms) with the same 
normalization as the detector signal and the nominal shapes of the distributions (red curves). 



distributions as is done with the two leptons. We therefore propose to use consistency cuts, 
as discussed in to purify the events. We assume that the endpoints of the distributions 
have already been estimated, but not necessarily very precisely. To plot the distribution 
of a given invariant mass among the set: rnq^i^^^, ?7iq;(jjigh) ^qlh we then cut away all 
events except those where one and only one of the two hardest jets, when combined with 
the leptons, gives invariant masses that lie below both of the endpoints of the distributions 
we are not plotting, i.e. that the jet we pick for an event in a given distribution is consistent 
with the other distributions and that there is no other such consistent jet, among the two 
hardest ones. Additionally we require that the two leptons have an invariant mass below 
the endpoint of the rriu distribution. As we shall see, these consistency cuts are very 
effective in leaving only events where we have picked the right jet to go with the leptons. 
However they will reduce the number of accepted events significantly, and thus increase 
statistical errors. 

The distributions of the invariant masses ?ti^;(1o„), m^;(high) a-^id niqii, after detector 



simulation, after the cuts of section |3^ and after cuts on 6-tagged jets and consistency 
cuts, are shown in figure ^ with error bars (black). For the consistency cuts we have used 
as endpoint values m^^^^^ = 320 GeV, m'^if^-^^^ = 395 GeV and m^f^ = 450 GeV.^ These 
lie ~ 20 GeV above the fit values of |12], and are clearly well outside the level of uncertainty 
expected for the measurement of the endpoints, showing that the consistency cuts will not 
critically depend on measuring the endpoints very accurately. Also in figure ^ (blue his- 
tograms) we show the parton level distributions after cuts, rescaled to the reduced number 
of events in the detector level distributions. Finally, we show the nominal distributions, as 
given from our formulae (red curves). 

From figure ^ we see that the distributions for mg^^ow); "T-(?/(high) aiid niqii in general 
remain unchanged after detector effects, but compared to figure ^ we have fewer events 
and thus larger errors. In the mg;(high) distribution we see a rounding of the peak due 
to additional smearing from energy measurements in the detector simulation. All three 



*For m™"^ we use 80 GeV which is ~ 3 GeV above the fit value. 
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Figure 10: Example of a "foot" and a "drop". 

distributions are also shiftet slightly to the left (lower invariant masses) of the nominal 
shape as a result of inperfect jet recalibration. For the niqii distribution we note that 
some events with an erroneously picked jet remain to the right of the nominal endpoint of 
the distribution. The number of these events can be reduced by tighter consistency cuts. 
Knowing the expected shape of the signal distribution, one could also model the shape of 
the background and subtract it. 

5. Feet 

As discussed in the mass values may be such that some invariant mass distributions 
exhibit "feet" or "drops" at the high ends. These can be hidden by a significant presence 
of background, taken to be smearing from detector effects or even assumed to be a width 
effect of the sparticles, making a precise determination of the kinematic maximum, and 
through these the SUSY masses, difficult and subject to systematic errors. 

There are two basic features in the distributions that can result in feet: The first is that 
the last function piece has a maximum value that is much lower than the global maximum, 
and so is not taken to be part of the distribution. This we will refer to as a "foot" . The 
second feature comes about when the last function piece does not fall gradually to zero, but 
ends in a discontinuous jump. This we will call a "drop". The two situations are illustrated 



in figure |10|. We will often refer collectively to "feet" and "drops" simply as "feet" . 

For both low value of the ratio between the maximum of the last function 

piece or the height of the drop, and the global maximum of the distribution will indicate 
the possible danger of mismeasurements due to feet. We denote this ratio r: 

height of "foot" or "drop" 

r = rr-^ : • (5-1) 

global maximum 

Without reference to a specific model and thus the size of the background compared to 
the signal, one can not give an exact value of r where this danger is real, but what can be 
done is to perform a model independent exploration of what relationships between mass 
parameters give low ratios. 

We will here discuss this for the invariant masses rn^2{high)j "^c2(iow) rricba^ in terms 
of the mass-squared ratios Ra, Rb and Rc introduced in ( p. 34 ). With experimental data 
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one could then use knowledge of the signal to background ratio, with measured values of 
the mass-squared ratios^, to look for such mismeasurements. Finally we will look at the 
situation for these invariant masses in the mi/2~"^o-pa'ra'meter planes around the Snowmass 
mSUGRA benchmark points SPSla, SPSlb, SPSS and SPS5 fll. 



5.1 mc2(high) 

We start our discussion of feet in the different mass distributions by an overview of the 
number of function pieces involved in each region. For the two-particle masses mc2(high) 
and ?n-c2(iow); these are given in table ||. 

For the ?n.c2(high) distribution the 

mass can be factored out of all 
the function pieces for any value of 
mc2(high)- Taking the ratio of two 
function pieces, these factors cancel, 
making the ratio independent of Rc- 
For the first two regions we can then 

look for a possible foot by finding the ratio r of the maximum of the last function piece to 
the global maximum in terms of Ra and Rb- The third function piece of the third region 
does not end at zero, thus it may have a dangerous drop, as described above. The ratio r 
of the function value at the endpoint to the global maximum can likewise be given in terms 



Distribution 


Region 1 


Region 2 


Region 3 


"^c2(high) 


3 


3 


3 


"^c2(low) 


1 


2 


2 



Table 1: Number of distinct functions for two-particle 
masses. 



of only Ra and Rb- The r-values of all three regions are plotted in figure 11. The results 
largely agree with those obtained for a particular set of mSUGRA model points in sec. 4.1 
of [l^, where the shape of the distribution was generated by a Monte Carlo decay routine. 

We find that for Region 1 the ratio is large for most values of Ra and Rb, so that the 
danger of a foot is small, unless we have a large background to signal ratio. On the other 
hand, in Region 2, there are large areas with a small value of r where we can have foot 
effects. One can show that, as a function of Ra and Rb, r is for ?nc2(high) continuous over 
the boundary between Regions 1 and 2. Approaching the border to Region 3 the ratio goes 
to zero. In Region 3, r again takes on large values, even more so than for Region 1. A 
minimum around r ~ 0.60 excludes the drop scenario from being a danger in this region. 

5-2 mc2(iow) 

For the ?Tic2(iow) distribution we can also factor out m/j in all regions, with the same factor 
for all function pieces, so that r again is independent oirriD- In Region 1 the only function 
piece does not end at zero, and we may thus encounter the drop situation. In Regions 2 
and 3 we have a potential foot from the second function piece. In figure O, we plot the 



ratio r of eq. (|5.lD . 

In Region 1 we find a strip along the border to Region 2 where we can have a dangerous 

drop in the distribution. Indeed SPS la lies in this region, which is somewhat surprising 

given that this feature was not noticed in |jl^ . This clearly shows the advantage of having 

analytical expressions for the shapes of the distributions. However, the drop here is so 

small, r ~ 0.046, that the resulting mismeasurement of the endpoint amounts to 0.1%, 

much less than the statistical error. For Regions 2 and 3 r is again continuous over the 
^The ratios Ra, Rb and Rc can to some extent be determined from the shapes of the distributions 
alone, and are thus not so susceptible to the effects of a mismeasurement of an endpoint. 
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Figure 11: Contour lines of r for all three regions of the mc2(high) (left) and mc2(iow) (right) 
distributions, plotted in the i?yii?B-plane. We also show the position of the points SPS la (a), 
SPS la (/?), (i), (ii) and (iii), referred to in the discussion of "feet" in 

border. The potential for a foot is however restricted to a small strip along the border to 
Region 1. 



Region 


Subreg. 1 


Subreg. 2 


Subreg. 3 


1 


3 


3 




2 


3 


3 




3 


3 






4 


4 


4 


4 



5.3 rricba 

The numbers of function pieces found in each region and subregion of the three-particle 
mass distribution are shown in table ^. 

In the three-particle case we can no 
longer remove Rc from the discussion 
of feet. However, other simplifications 
arise. One can show that none of the 
four regions of the rricba distribution has 
a drop as the last function piece goes to 
zero at the kinematical endpoint for all 
regions and subregions. What remains 
is the possibility of a foot. In figure |l^, 
the ratio r is shown for two values of Re- 
in the entirety of Regions 1, 2 and 3 we have r = 1 with no danger of feet because 
the global maximum is in the last function piece. Only in Region 4 can a foot occur. The 
dangerous regions are near the border between Region 4 and either Region 1 or Region 2, 
which have low values of r and potentially a misleading foot. Comparing the plots for 
Rc = ^ and Rc = | we see how the value of r changes with Rc inside Region 4, and in 
particular observe that the points where Regions 1 and 2 meet and where Regions 1 and 
3 meet also depend on Rc- However, for rricba, the final function piece often has a steep 
negative slope after the maximum, which reduces the possible negative effect of a foot. 



Table 2: Number of distinct functions for the three- 
particle mass rricba- 
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Figure 12: Value of r in all four regions of the rricba distribution, plotted in the i?^i?B-plane for 
Rc^ ^ (left) and Rc ^ ^ (right). 



5.4 SPS benchmark points 



In figures 13 and 14 we show the value of r in the mi/2~™'0-pla'nes around the Snowmass 
mSUGRA points SPSla, SPSlb, SPSS and SPS5, for the mc2(high) a-iid mc2(iow) distribu- 
tions respectivly. In the top left panel we have tan /? = 10 and = —^o and in the upper 
right panel, tan/3 = 30 and = 0. The lower left panel has tan/3 = 10 and Aq = and 
the lower right panel, tan /? = 5 and Aq = —1000 GeV. We have only considered parameter 
values where we have decay chains of the type given in eq. (U), with an on-shell, right- 
handed slepton. For the analogous decay chain via left-handed sleptons, there is no danger 
of feet for values of mi/2 up to 1 TeV. In the gray area we have m^o < ^i^, thus the 
decay is only possible via a virtual slepton. The white area, marked TF, is theoretically 
forbidden and the light brown area has a charged LSP. 

In the mSUGRA planes we see that the dangerous area for the 7Tic2(high) distribution 
lies in the narrow Region 2, and for the mc2(iow) distribution in a fairly narrow corridor 
along the border between Regions 1 and 2. The WMAP-consistent bulk region at low 
aud mo in the top left panels of figures and U, around the point SPSla, lies in 
these dangerous areas. This can be compared to the SPS3 parameter line running along 
the WMAP-consistent stau coannihilation region in the bottom left panel, and we see that 
there is little risk of problems with feet in the stau coannihilation region for values of mx/2 
that are still viable when we consider LEP limits on SUSY masses and on the lightest Higgs. 
However, this conclusion is only valid for relatively low values of tan/3. In the top right 
panel we have tan/3 = 30, and indeed the dangerous area crosses the stau coannihilation 
region running along the border to the region with a stau LSP. The lower right panel 
shows that the WMAP-consistent stop coannihilation region, to the left of SPSS, along the 
theoretically forbidden region, is also dangerous in this respect. 

We do not show results for nicba in these parameter planes as we find only points in 
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Figure 13: Value of r for the mf,2(high) invariant mass distribution in the TO]^/2~"^o-planes around 
the Snowmass benchmark points SPSla, SPSlb, SPSS and SPSS. 

Region 1 of the nicba distribution in the mSUGRA planes we consider here, for values of 
mi/2 < 3 TeV. As discussed above, Region 1 in this distribution contains no dangerous 
"feet". 



6. Summary 

We have derived analytical expressions for the invariant mass distributions of massless 
SM endproducts c, b, a, in cascade decays of the form D Cc Bbc — >■ Aabc. Our 
main results are valid for the decays of spin-0 particles, D, C, B, or equivalently, for a 
sum over all combinations of charge and chirality in the final states. In a hadron collider 
environment, the difficulty of determining the charge of quarks from the jets they instigate. 



- 31 - 




Figure 14: Value of r for the 'm^2(iow) invariant mass distribution in the rrii /2-Tni^-"p\a.Ties around 
the Snowmass benchmark points SPSla, SPSlb, SPSS and SPSS. 



makes this a reasonable simplification. We have discussed how different spin configurations 
can easily be included in the distributions, and in the appendix show the corresponding 
distributions for the specific SUSY decay chain given in (|l.lD . 

The effects of cuts, particle widths and final state radiation on the shape of the distri- 
butions have been investigated with PYTHIA in a particular SUSY scenario, the Snowmass 
benchmark point SPSla. We find that while a set of cuts used to remove SM background 
have the potential of distorting the expected shape we can identify which cuts are re- 
sponsible for this, and in which regions of invariant mass this distortion takes place, by 
looking at the distribution of the cut parameters, and without resorting to Monte Carlo 
truth information. The effect of particle widths can be compensated for by a parametrized 
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smearing of the distibution. Final state radiation introduces a systematic loss of energy, 
thus shifting the invariant mass distributions towards lower values. 

We have also studied the effects of a generic LHC detector, in the same SUSY model, 
through the use of the fast detector simulation AcerDET. The analytic expressions for the 
invariant mass distributions is seen to survive the smearing effects of the detector, and we 
demonstrate a method for handling the combinatorial problem of picking the correct jet 
that belongs to the decay under investigation, from amongst the many candidates. This 
consistency cut method is shown to be very effective in removing combinatorial background, 
but results in a significant reduction in the number of events. The question of whether it 
is optimal to use cuts to reduce this background, or whether to try to model and subtract 
it, is still open. A small, systematic shift of the distributions towards lower invariant 
masses is observed. This indicates that the reconstruction of jets and the jet recalibration 
routine in the detector simulation, is insufficient for jets from this decay chain, and imply 
the unsurprising conclusion that understanding the jet energy scale will be essential in 
reducing systematical uncertainties in parameter determination from the shapes of the 
invariant mass distributions. 

We also demonstrate an application of the analytical expressions in finding dangerous 
"feet" in the distributions, that could lead to mismeasurements of distribution endpoints, 
and in turn masses. In scans over mSUGRA parameter space we find that this danger exists 
for several WMAP-consistent regions: the bulk region and the stau- and stop-coannihilation 
regions. 
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A. Including spin in a SUSY scenario 

One can add spin effects in the decay chain by multiplying the integrand with a suitable 
angular dependent function fuiu) or fv{v), depending on the spin configuration in question. 
We will here give the results of including spin in the SUSY decay chain (1.1), as discussed 



in the Introduction. In this decay chain we have A = Xi, B = I, C = X2 D = qi. The 
spin-^ of X2 '^ill yield an extra factor of either 2u or 2(1 — u) depending on the different 
combinations of chirality and charge in the final state for the quark and the "near" lepton, 
as given in table ^. Only the decays starting from a left handed squark are listed fully, those 
for the right handed squark follow from a simple interchange of left and right handedness^. 
The "far" lepton does not contribute any spin-dependent factor, since the LSP, associated 
with that fermion line, is not observed. 



'Note that the handedness of sparticles refer to the chirahty of the SM particles they couple to. 
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R 
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L 
L 
L 
L 
R 
R 
R 
R 
L 

R 



chirality of If 



L 
L 
L 
L 
R 
R 
R 
R 
L 

R 



Factor 



2u 
2(1 - u) 

2u 
2(1 - u) 
2(1 - u) 

2u 
2(1 - u) 

2u 
2(1 - u) 

2u 



Table 3: Spin factors. The chirality L/R of Z+ (q) denotes that it is the antiparticle of a left-/ 
right-handed (q). 



Except for the first case, discussed in section A.l, we shall not give explicit formulas 



for both chirality cases. The two cases are related as follows: 



1 dV 



r dm? 



2 dVo 1 or 



2(1-m) 



Tq dm? r dm^ 



(A.l) 



2u 



The integrated decay widths F are for both chirality configurations equal to Fq, since 
the spin correlations only introduce forward-backward asymmetries (in the X2 '^^st frame), 
which integrate to zero. 



A.l mr 



We begin by including spin in the mca distribution by multiplying the RHS of eq. ( 2.14 ) 
with fu{u). We find for the /«(«) = 2u case that 
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(A.2) 
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(A.3) 

This confirms the result of [20|. As required, the average of the two again gives eq. ( |2.2C| ). 
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A.2 mc2(high) 

For Region 1 of the mc2{high) distribution, multiplying ( |2.26| ) by fu{u) = 2u gives 
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(A.6) 
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A.3 mc2(iow) 

For the mj,2(iow) distribution, in Region 1, and with = 2m, we get 
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for < mc2(iow) < m^^. For Regions 2 and 3 the distribution is for fu{u) = 2u given by 
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"^c2(low) 



+ 



(^max-)2 (m^axp_^^2 



(^max)4 

2 2 
_^ '^c2(low) _^ ""^c2(low) ^ 



^max 



for 



mc 



c2(low) 
< "ic2{low) < "i:;^2(eq) 



(Ai 



A.4 mcfea 



In the mcba distribution the spin factor enters into the integral L{ai, a2) of eq. ( 2.71| ). We 
define a new integral: 



ra2 

M{ai,a2) = / 

J ai 



y 



zdy 



y"^ + 2ymB + m^ — mj-^ 
y2 + 2ymB + m^ — m|) 



mBL{ai,a2). 



(A.9) 
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For all four Regions of the rricba distribution the effect of the neutralino spin can be added 
by making the following simple substitution in the distributions of eqs. ( p.89[ )-( ^.96 ), for 
the fu{u) = 2u case: 
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